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Abstract 


The Geogiaphic Information Systems ( G I.S ) are a powerful set of tools for collect- 
ing, storing, retrieving at will, transforming and displaying spatial data fioin the real 
world for a particular set of purposes Presently it is increasingly used by the hydrol- 
ogist for spatial data analysis and decision making GRASS (Geographical Resources 
Analysis Support System) is a raster based G I S. developed by USACERL and is 
available for study The additional components of this software are implemented at 
SPARC LX workstation at Water Resources Engg and Management centre, Depart- 
merit of Crvjl Eiigineeniig, 1 1 T Kanpur. Then GRASS was used for the spatial data 
base development of Ganga River Basin between Kanpur and Allahabad and for thi* 
ground water balance analysis of this region The maps are prepared according to 
the ipqujicments, from the maps collected from U.P. Irrigation Department. These 
maps are scanned and then digitized using the on seteen digitization capabilities of 
GRASS {v. digit). Then the raster maps of the study region. Rain gauging stations. 
Districts m the study region, Ground water measurement stations, Irrigation Blocks, 
Channel Networks, Main rivers etc , are prepared using the program v.to.rasL 

The normal rainfall values for monsoon and nonmonsoon season aie calculated 
at each rain gauge station and are implemented in tire data base. Further the.se 
rainfall values are used to interpolate the rainfall distribution in the study legion. 
this interpolation is done using the utility r.surf.idw. Similarly the pre monsoon and 
post monsoon ground water table values at different monitoring stations are collected 
and are used to interpolate their distributions over the entire study region 

From the topographic map of the study area the topographic elevations at dilferent 
points afe estimated using the program r.surf.idw. The ground level contours are 
drawn from this interpolated surface. 


II 



(yontouis of rainfall and giouud vvatet table levels for both the seasons and the 
topogiapluc contours aie extracted from their interpolated surfaces. The depths to 
giound water table at diffeint iingation blocks are calculated. 

The average rainfall values in different districts, irrigation blocks and nr the entue 
study area is calculated for the two seasons Diffeient components of giound water 
balance aie calculated and irsing a lumped model approach an appioximate ground 
water balance study has been earned out for the study region. 
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Chapter 1 


Introduction 


1.1 General 

Tiaditionally decision making process in Water resources management and planning 
has been assisted by a multipbcity of decision siippoit tools, including mathematical 
models for simulation and optimization, statistical routines for analysis and interpie- 
tatioii of data and database management systems Recently with the advent of ever 
increasing computer power, Geographic Information Systems(GIS), have been added 
to the body of the tools that are intended to support the decision processes thiough 
their many steps. 

Today large Geogiaphical Information Systems aie being developed, opeiated , 
continuously updated and furthei refined by various environmental and water agen- 
cies. This emerging new technology was developed at first mainly fiom multiobjective 
representing, structuring and mapping of multivariate spatial data The GIS potential 
is now being realised not only m inventory application but also in studies which apply 
simulation models. Tailored to a special problem region of interest and confined to 
the set of problems relevent data and models, a GIS can be applied both as a data 
platform for simulation model as well as an interactive user environment for daily 
routine water management in smaller agencies. In many cases, however, GlSs are 
laid out to serve as multipurpose, nearly ’’catch all” data management .systems with 
numerous possible applications, ranging from hydrological and environmental studies, 
through land use and economical investigations up to support of urban development 
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Figme I 1: Database Management in G\S (Valenzuela) 
and transpoit planning. 

Leanfear's{[%9) editorial ui Water Resources Bulletin concludes that, "The im- 
pact of GIvS may be as significant as the introduction of FORTRAN, as scientists 
learn to use and manipulate spatial data to contiibute to new level of understanding 
of environinontal issues”. Although not enough time has been elapsed to accurately 
assess this claim, the potential application of GIS to Water Resouices is indeed great 

GIS are used to assist decision makers by indicating various alternatives in devel- 
opment and conservation planning and modeling the potential outcomes of series of 
sceneries as shown in Figure 1.1. 

Computeiized mapping and spatial analysis has been developed simultaneously in 
sevcial related fields such as, 

• Cadastral and topographic mapping, 

• Thematic cartography. 

• Civil Engineering 
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• Geogiaphy 

• Hydrology etc. . 

1.1.1 Typical GIS applications: 

• automatic cartogiaphy 

• subdivision design(cut/ftll. street and paicel layout) 

• cadastrel mapping 

• highway mapping 

• utility and facility mapping and management. 

• geodesic mapping 

• event mapping 

• census and related statistical mapping 

• management of well log data 

• land use planning and management 

• maiketing and retailing 

• vehicle routing and scheduling 

• uibaii and regional planning 

• route selection of highways, pipelines and similar data 

• water resources planning and management 

• suiveying and eiigineeiing 

• mineral exploration. 

Current types of G.I.S.s ; 

Current types of CIS are grouped into; 
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• engineeiing mapping systems (CAD/CAM) for applications in photogrammetry, 
topographic base maps and road engmeeiing. 

• piopeity of paicel infoimation systems 

• generalized thematic and statistical mapping systems toi natural lesources man- 
agement, forest inventories, vegetation, geology, soils, census mapping 

• geogiaphic base file systems associated with street networks 

• image piocessmg systems 


1.2 Definitions of G.I.S. 

Theie have been a number of attempts to define a Geographic Information Sys- 
lems(GIS), and two mam themes can be distinguished, technological and problem 
solving The technological definition typified by the followings, concentrates on 
the computei related aspects of the field, and has dominated the liteiature till to 
date (^Carre 11 , 1984, Marble, 1984,Pci’f'ker, 1988) 

Definition 1: 

Geographical Information Systems are a powerful set of tools for collecting, stor- 
ing, tiansformiiig and letrieving at will and displaying spatial data from real world for 
a paiticular set of purposes (Burrough, 1986, p.6) 

More recent definitions have concentrated on problem solving aspect of GIS. 
Tims Goodchild (1985) gives the following definition. 

Definition 2: 

A GIS is best defined as a system which uses a spatial database to provide answers 
to queries of a geographical nature. 

This emphasizes the analytical nature of the system. An emphasis on problem 
solving nature has been presented by Cou)en(l987) 
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Definition 3; 

A GIS IS a decision suppoit system involving the integration of spatially referenced 
data m a pioblem solving environment (Cowen 1988 p 1554) 


1.3 Inter Relationships Among Cartography, Re- 
mote sensing and GIS 

The model illustrated m the Pig, 1 2 is presented as the most leahstic representation 
of the intei action among the three fields as currently practiced No field is placed in 
a position of dominance one over the other, oi indeed in isolation, but there is intei- 
action in all possible combinations of the tliiee Data acquisition and data analysis 
aie the emphasis m remote sensing. With m lemote sensing however , considera- 
tion must be given to the destination of the iiifoimation gathered, which may either 
be used as input to a GIS or be repiesented m a graphic form. Furthermore, the 
computer system used in modern digital image analysis have much m common with 
GIS processing, such that many vendois have, without changing the data structures, 
implemented GIS function in their image processing system Clearly, then, remote 
sensing does impinge on the domains that might otherwise be recognized as GIS and 
cartography, respectively 

The analysis of geographic information to support decisions, which is the main con- 
cern of practitioneis of GIS, is dependent on the way in which the data are gathered. 
When data are derived from remote sensing, quality and organization are dependent 
on the methodology in that field. Similarly while GIS professionals aie not primarily 
conceined with the quality of the graphics that may be derived from the informa- 
tion, they should be cognizant of the implications of the data manipulations upon 
the message presented in a resultant map. The cartographer may map information 
which is a direct product of remote sensing or which have been processed by a GIS. 
In this model, the emphasis within cartography is the effective presentation of the 
information as a map, and includes data analysis and manipulations to facilitate that 
presentation. 

The relationship of remote sensing and GIS with cartography in this model is not 
unlike the traditional view of data acquisition and data manipulation as precursor 
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Figure 1.2* Model of three way interaction of Remote sensing, Caitography and 
G.I S(Burrough) 

of the map The three way interaction model, however, depicts a relationship that 
recognizes laigei siibsideiies of caitogiaphy Remote Sensing and GIS have a body of 
knowledge and lequire a level of expertise such that it is no longer realistic to envisage 
a cartographer as being expert in these aieas as well Thus the lelationship of tlie 
fields requites not only that the definitions are reasonable but also they reflect the 
nature of the interactions. A set of revised definitions are suggested 

Cartography; 

Cartography is the field which is involved with the graphic communication of 
spatial relationships and distributions, and includes the analysis and manipulation of 
geographic data to enhance representation. 

Remote Sensing; 

Remote sensing is the capture and interpretation of data from regions of tlie 
electromagnetic spectrum through the use of noncontact instruments together with 
analysis and manipulation to facilitate interpretation. 

Geographical information systems; 

GIS is defined as the management analysis and manipulation of spatially refer- 
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oiiced iiitoimation in a problem solving synthesis By observing the above definitions, 
Analysis and manipulation appears to be forming the common theme 

In caitogiaphy analysis and manipulation involves in the piocesses of projection 
change and line generalization and finding the chloiopleth map intervals. In remote 
.sensing cluster analysis of reflectence forms the analysis and manipulation pait And 
in GIS oveilay analysis and Zone buffeiing forms the iinpoitant phase of analysis and 
manipulation 


1.4 Components of G.I.S. 

Geographical Infoimation Systems have three important components, computer hard- 
ware, set of application softwaie modules and a proper organizational context. These 
three components need to be in balance if the system is to function satisfactorily. 

Computer hardware components of a GIS are presented in Fig.l 3 The computer 
02 the central processing iinit(C.P.U.) is linked to a disk drive storage unit, which 
provides space for stoiing data and ptograms A digitizer or other device is used to 
conveit data from maps and documents into digital forrei and send them to computer 
A plotter or other kind of display device is used to represent the result of data 
processing, and a tape drive is used for storing data or programs on magnetic tapes, 
or for communicating with other systems The user controls the computer and its 
peripherals via a visual display unit(V D U), otherwise kirown as a terminal. The 
user’s teiminal might itself be a micro computer, or it might incoipoiate special 
hardware to allow maps to be displayed quickly. 


1.5 G.I.S. software modules 

The software package tor a GIS consists of five basic technical modules (Fig. 1.4). 
These basic modules are subsystem for : 

• Data input and verification 

• Data storage and database management. 
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Figure 1.3 Hardware components in a GlS(Burrough) 



Figure 1.4: Main software components in a GIH( Burro ugh) 
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• Delta output and piesentation 

• Data tiansformation 

• Intel action with the user 

Data input(Fig 1 6) coveis all the aspects of transforming data captured in the 
foim of existing map, field observationa and sensois( including aeiial photogiapliy, 
satellites and recording instruments) into a compatible digital form A wide range 
of computei tools aie available for this puipose, including the interactive terminal 
01 Visual Display Unit (V D U) The digitizer lists the data in a text file, scanneis 
(in satellites or aero planes direct lecordiug of the data or foi converting maps to 
photographic images) and devices necessary for lecording data already written on 
magnetic media such as tapes, diunis and disks. 

Data stoiage and database managemertt(Fig 1 5) concerns the way in which the 
data about the position, linkage( topology) and attributes of geographical elements 
(points, lines, and areas lepiesenting objects on earth’s suiface) ate stiuctured and oi 
ganised, both with respect to the way they must be handled in the computer and 
how they aie perceived by users of the system The program used to organise the 
database is known as database management system (DBMS) 

Data output and piesentation(Fig 1 6) concerns the ways the data are displayed 
and the results of the analysis are reported to the user. Data may be presented as 
maps, tables and figures in a variety of ways ranging from the ephemeral image on a 
Cathode Ray Tube(CRT) through hard copy output drawn on printei or plotter to 
information recorded on magnetic media in digital foim. 

Data transformation{Fig.l.7) embraces two classes of operation, namely (a) trans- 
formation needed to remove errors from the data or to bring them up to date or to 
match them to other data sets and (b) the large array of analysis methods that can 
be applied to the data in order to achieve answers to the questions asked to the GIS 
Tiansformation can operate on spatial and non spatial aspects of the data, either 
sepal ately or in combination. Many of these transformations, such as those associ- 
ated with scale-changing, fitting data to new projections, logical retrieval of data and 
calculation of area and perimeters are of such general nature that one should expect 
to find them in every kind of GIS in one form or another. 
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Figure 1.7. Data Transformation ^Swrrou^/iJ 



Figure 1.8: Organizational aspect of a GIS (Burrouffh) 

The five technical sub-systems of GIS govern the way in which geographical mfoi- 
mation can be processed but they do not of themselves guarantee that any paiticular 
GIS needs to be placed in an appropriate organizational context as in Fig.l 8. 


1.6 Capabilities of G.I.S. The four Ms: 

Ones understanding of the planet has always been limited by the lack of information, 
as well as lack of wisdom and knowledge. For things too small to see, one has devel- 
oped microscopes that can image down the molecular level.At the other end of the 
continuum for things that are (in very real sense) too large to see, one has geostatioii- 
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Figure 1 9. The four M’s m a G1S( Burro ugh) 

ary satellites that can take an image of an entire hemispheie Geogiaphic Information 
Systems are a means of integrating spatial data acquired at different scales and times 
and in different formats. Basically urban planners, scientists, resources managers and 
others who use geographical infoimatioii, work in several main areas. They observe 
and measure environmental parameters They develop maps which portray charac- 
teristics of the earth They monitor changes in human surroundings in space and 
time In addition they model alternatives of actions and processes operating in the 
environment These, then, are the four Ms: Measurement, Monitoring, Modeling and 
Mapping{Fig,l.9). These key activities can be enhanced through the use of informa- 
tion system technologies and in particular , through the use of a GIS. 


1.7 Different types of G.I.S. 

As mentioned above a GIS can be described as an information processing computer 
system that handles and analyses internally referenced spatial and non spatial data. 
Primary GIS operations are: e.g, thematic mapping, area and distance measurement, 
and spatial overlay. Sequences of primary GIS operations can be linked by the so 
called compound GIS operations for spatial analysis purposes. Figl.lO illustrates 
some GIS features and outputs 
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Figuie 1.10: Examples of GIS data representation and object manipulationfflurrouij/i^ 
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Theie aie two ctifFeieut types of GIS • 

• Rastei based GIS 

• Vector based GIS 

In a raster GIS data are stored and manipulated as picture ele]Tients(pixels) of reg- 
ular two-dimensional matrices (digital images, bitmaps), stacks of multiple map layers 
are used to provide multivariate data representation of a region These GIS versions 
strongly support the representation of pictorial information such as multispectral dig- 
ital images from remote sensing systeins(Fig 1 10) as well as the implementation of 
digital image processing operations(e.g density slicing, spatial filtering, multispectral 
analysis). Therefore, many commercially available raster GIS are offered with options 
to incorporate digital image processing modules 

A vector GIS generally features relational database management, object oriented 
data representation as vector graphics(nodes,Iine segments and polygons Fig 1.10 and 
numerous object oriented search and analysis operations (e g Boolean query,geometry 
analysis and measurement functions, topology analysis). In addition to area covering 
spatial information, vector based GIS can handle non spatial textural, temporal and 
tabular information The linkage of GIS and a simulation model is well supported 
by object oriented data representatioii(Fig 1 lOd Arnold and Orlob, 1988) Advanced 
commercially available vector GIS packages( ARC/INFO) consist of four mam com- 
ponents; the dialog module(interactive user interfaces), the relational database man- 
agement system, the data analysis and manipulation module and the data display 
module. A comparison of both raster and vector methods are given in Tab. 1.1. 


1.8 Objective of the study 

The objective of the study is to implement additional components of already imple- 
mented (partly) GRASS4.1, a GIS software GRASS4.1 in the Sun-SPARC LX work- 
station of Department of Civil Engineering, I.LT Kanpur, to test the XGRASS4.1, an 
interactive version of GRASS4,1(GUI) in the same system, to illustrate it’s capabili- 
ties by creating the hydrologic data base for the Ganga River Basin between Kanpur 
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RASTER 

VECTOR 

ADVANTAGES 

1) Simple Structure 

2) Easy Oveilaying 

3) efficient for high 
spatial variability 

1) inore ’’compact” structure 

2) efficient topology encoding 

3) better suited for graphics 
{eg hand-drawn maps) 

DISADVANTAGES 

1) Less compact(if not compressed) 

2) topology more difficult 

3) Jaggedness in output graphics is 
seen unless files are veiy laige 

1) More complex structure 

2) overlaying more difficult 

3) Inefficient foi high 
spatial variability 

4) not effective for ehancement 
of digital images 


Table 1 1 Comparision between Raster and Vector Methods 

and Allahabad and utilizing this data base for the ground water balance study of the 
area in both the monsoon and nonmonsoon seasons. 


1.9 Scope of the Study 

This study is not designed to develop procedures and software, but to utilize the 
existing software efficiently. The purpose of the study is to address the development 
of a GIS focusing the role of automated data processing in improving the effectiveness 
and efficiency of data management and data analysis using the available data for the 
study region 


1,10 Significance of the Study 

In the past, relevant data have been available from a variety of map sources. Though 
visually accessible, these map sources do not permit easy comparisons with site data 







[ntrodiLction 


16 


and most impoitantly, they do not permit the essential considerations of two oi more 
data sources. A variety of computeiized systems aie now available that have these 
multivariate capabilities to deal with complex spatial data These GIS are conceptual 
equivalents of other computerized database management systems in that they peimit 
a generalized automated stiuctiue for the entiy, storage and retrieval of geogiaphic 
data Though similar to other database management systems, they have the added 
complications presented by the spatial component of the data and its substantial size. 
A final difFeience is that GIS requites substantial commitments to data gathering m 
oidei to obtain information on the geographical chaiacteiistics of a particular area 

Ceilainly a technology is not a panacea for all ills. It is hoped that the effective 
and iiitegiated database management GIS technology, lemote sensed imagery and 
appropiiate predicative modeling can enhance management of resources 


1.11 Organization of the study 

The study is repoitecl in following sequences- 

• Introduction of GIS, objective, scope, significance and organisation of the study 
aie discussed in Chapter 1. 

• Desciiption of the GIf3(GRASS) is discussed m Chapter 2. 

• Development of a spatial database is discussed in Chapter 3 

• GIS features of GRASS, extraction of rainfall contours (for both monsoon and 
non monsoon seasons) as well as ground water table contours( both for pre- 
monsoon and post-monsoon) , average rainfall using Theissen polygons, and 
database creation for different characteristics of irrigated aieas in the study 
region are discussed in Chapter 4. 

• Ground water balance for the two seasons is presented in Chapter 5. 

• Summary, conclusions and suggestions for fuither study are presented in (Chap- 
ter 6. 


Chapter 2 


Geographical Resources Analysis 

Support System 


2.1 Background 

The Geogiaphical Resoiiices Analysis Siippoit iSystem (GRASS) is <a niiciocomputei 
based image pi Dressing and Geographic Infoi ination Sysfcem(G IS ) used to develop, 
manipulate, analyze and display geogiaphical data sets Initially developed by u'- 
seaicheis at USACERL’s Enviioniuental Division foi envitonmental planners at mili- 
tary installations, GRASS is now used by a vaiiety of public and puvate agencies ainl 
individuals to access environmental impacts, evaluate site suitability, detect change, 
manage resouices and model the effects of environmental phenomenon across a land- 
scape. 


2.2 Definition of GRASS 

Ditfeient G I.S have dilfeicnt ohjectivc',, abilities, opeiatmg e'liviionments and siojn' 
Mo single system does eveiy thing. (tlt/VSS ( an lie Itiiedv delmed by noting ili.it il 
is. 

A Geographical Information System(G.I.S.) 
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It IS not a word processor, a statistical system, a data base management system, 
or ail automated map ping- facility management system It is possible to interface 
GRASS to such systems however 

A Program that accommodates raster and vector data format - 

Both laster and vector data foimats aie accommodated (to different extents) ui 
GRASS Vectoi data are used in GRASS for digitizing and the graphic overlay of 
data AH data analysis is done on raster data, and GRASS is therefore commonly 
leteied to as a raster based system GRASS contains coiiveision programs to translate 
vector data into a laster foimat and vice-veisa. 

A Set of tools 

When arianged with other systems, GRASS can be employed to successfully man- 
age land lesoutces 

A data input-output system 

GRASS IS capable of’ data digitizing, data lead-in and data lead-out conveisions, 
image piocessing, data analysis and data presentation GRASS continues to glow, 
and to piovide moie sophisticated utilities (loopeicdive effoits coutimie to be need 
ed to cooidmate the use ot GRASS with the system, technologies, ijersoniiel, and 
regutatious alieady involved in the planning piocess 


2.3 Data and Capabilities 

The way in which computei data is oiganized defines the potential for its use Geo- 
giapliic data can be stored in two basic formats • uaster(grid cell), and vectoi (arc- 
notle) Image analysis systems typically use raster data, while computer aided design 
(CAD) packages typically use vector data. GRASS uses both data formats. 

Maps winch (ontatn ilistimt liiieai teatmes siuli as kmiIs <uu1 stioanis cU m 
distiiut leal teatuii's siidi <is roniiLiy houiidaiies, liaming aieas, sods, polygons, I’U 
aie input to GRASS in vector foimat. These featiues are defined and stoied a.s a 
sene.s of two dimensional cooiilinate paiis(polnls),this vector data is then i unverted 
into rastei foimat, since analysis ptogtams in GRASS work on lastei data. 
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GRASS uses the raster format for the processing of image and geographic data 
In a rastei format the landscape is divided like a checker board into a regular rectan- 
gular parcels of land attached to each parcel are identifying attributes specifying(for 
example) the particular paicels sod type, land cover, land use, vegetation, geology, 
slope, elevatioii,etc 

GRASS, at release of version 4 1, consisted of nearly 250 different computer pro- 
grams which the user can run directly through keyboard commands oi indirectly 
through menus and othei programs. A detailed list of commands are provided in 
Appendix These capabilities can be placed in following categories. 

• Geographic An<alysis 

• Image Processing 

• Map Display 

• Data input 


2.3.1 Geographic Analysis 

Grass provides several capabilities for map analysis and oveilay These include 

• Analysis of site data 

• Analysis of vector data 

• analysis of raster data 

• 2-D and 3-D manipulations 

• Single cell map reporting functions 

• Multi-Cell coincidence tabulations and comparisons 

• Neighborhood Analysis and filteis 

• legional analysis 

• Proxiimty analysis 
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• Boolean overlay functions 

• Ten am analysis 

It is tins capability of bringing data derived from satellite imagery, paper maps, 
and other computers(e.g. elevation data) together as input to some land use questions 
that peihaps best defines a G I.S. Maps originally at different scales and resolutions 
can all provide data which can be considered in some analysis 


2.3.2 Image Processing 

Aerial images are important data sources for a G I S. They come in the form of 
satellite images and high altitude photography 

Raw images contain tiemendous amount of information GRASS image process- 
ing tools provide two primaiy functions necessary to piepare data foi inclusion m a 
geographic information system. 

2.3.3 Map display 

Geographic and image processing generally requite that the operator have a significant 
amount of training to ensure useful output. Hence, these tools are not often of interest 
of the GRASS novice. The most attractive and useful tools for immediate use are 
image display capabilities. One set(the d tools)allows the user to manipulate the 
display of data on the computer monitor, while the other set (the p tools) enables to 
generate landscape images and maps on paper. 

2.3.4 Data input 

Perhaps the most impoitant components in a G IS. are those that aie used to captiiio 
the data GRASS data is derived from paper maps, satellites and other computers 
A powinfiil array of programs for reading magne^Lic tapes, manipulating raw data into 
a form usable with GRASS, and extracting information, allows for the production of 
quality data. The GRASS end users will likely become very familiar with dlgili'dng 
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utilities The process of entering mapped data into computer must be understood 
at every GRASS installation It piovides the user who depends on the data with a 
means of keeping the data current 

GRASS also provides an array of programs which will allow the user to read data 
from other sources Common sources include the Digital Line Graph(DLG) and Dig- 
ital Elevation Model(DEM) provided by U.S. Geological Survey and the Digital Ter- 
rain Elevation Data(DTED) provided by Defense Mapping Agency(DMA). Because 
data development is the most expensive component of establishing any G.I S., new 
programs are continuously under development to read and translate data originally 
created for other systems into a form usable by GRASS, 


2.4 XGRASS4.1 

XGRASS4 I IS a Graphical User Interface for GRASS4 1 and it offers the user, a 
tacility with which he can access the entire GRASS4 1 interactively. This software 
also has a GRASS shell through which one can use the grass commands also 


2.5 Portability of GRASS 

GRASS development has been accomplished on a variety of UNIX macliines writ- 
ten primarily in C language. It also makes strong use of UNIX system commands. 
UNIX is available on PCs, workstations, minis, mainframe and supercomputers, mak- 
ing GRASS relatively portable. 


2.6 Garden Concept 

Idic new level in the GRASS hierarchy has been dubbed the '’GARDEN” concept 
.Just one can put together many salads using items in vegetable garden, it is possible 
(o pick and choose from a huge assortment of computer programs m the analysis ol 
data GRASS provides only some of the items to choose from, the others are the 
standard programs provided with the operating system and may also include othcM 
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Figure 2 1. Garden concept in (Wasterlet) 


piopiietary and non-proprietary software acquired to help the user. The computei 
gaulen provides limitless opportunities for combining the computer fiuits While 
one mix GRASS with a Database management system (DBMS) and a report wiitei 
to assist enviionmental manageis to keep tiacks of oil spilling sites, another would 
mix GRASS with statistics package, a voice iccognition system, a natuial language 
inteiface and a DBMS to allow voice queiy of a forest management system(Fig 2.1) At 
gaiden level GRASS can remain a solid package of G I S. programs which provides the 
hoisepower to evei changing series of iisei interface and agrowing airay of computer 


programs 





Chapter 3 


Development of a Spatial Data 

base 


3.1 General 

G I S database development using GRAS’lS primarily m volves planning, building and 
maintenance of geographic data sets It can bioken clown into four stages. Planning, 
Data collection and map preparation, Installation of data into GRASS, and augment- 
ing and altering the existing data set. Both GRASS4 1 and XGRASS4.1 were used 
in the following steps for the study 


3.2 Planning 

A database is a collection of map layers on various themes for the same geographic 
location. A judicious decision is to be made regarding the selection of the layers of the 
database. This is followed by the spatial resolution which is needed for the map layer. 
A determination is then made as to which ot the debited map layers <ire available in 
digital or nondigital foimats at desired level of detail Finally, foi each map layet, the 
number of category types, (what level of detail will be represented) is decided upon 
All of these decisions will depend on the applications for which the map layers aie to 
be used. 
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Map layer prepatation and input to G I S. can be quite costly These naap layers 
necessaiy for G.I S. applications should be selected on a priority basis before the data 
ace acquired (while some map layers like basin boundary, hydrography and soils are 
fairly standard, several additional map layers can be derived from one or more original 
layers). Caieful planning coupled with technical know how and good understanding 
of the concepts, however decrease the data development cost to a significant extent 


3.2.1 Map scale and map layer resolution : 

Raster map layer data are stored in a grid, made up of tiny rectangles(cells) each 
cell in a grid represents an area on the earth which they represent, is represented by a 
resolution choosen for that map layer A map layer which has North-South and East- 
West lesolution of 50m, for example is a 50m*50m square on earth. Since larger scale 
map contains more information for a given place on the earth than does a smaller 
scale map, it is appropiiate to choose a coiiesponding cell size which lepresents a 
smaller piece ot earth{ in greater detail) for a larger scale map, and vice versa The 
lesolution used is leported in Tab 3.1. 

Resolution is expressed in a counter-intuitive way High resolution denotes greatei 
detail, but is expressed by small numbeis. Conversely, low resolution denotes low level 
of detail, but is expressed by large numbers. For example , 10m resolution data is of 
higher resolution than 100m resolution data.lOm resolution data is capable of showing 
100 times as much detail as is shown by lOOm resolution data. 


3.2.2 Cartography and digitization: 

The presentation of map for digitizing, and digitizing itself are very time consuming 
procedures that require careful attention to details, Utmost care need to be taken 
for geneiating accurate maps Gieat deal of tare and very .sinceie elf or t when put 
to digitizing pioces.s ran bring out accurate maps which would be further iisetl toi 
analysis. Besides having extensive hands on computer experience, one should also be 
able to read data typos, georeference satellite imageries, and exchange data between 
systems. 

Cai'logiaphy and digitizing task mdudes digitizing new or updated map informa- 



Development of a spatial Data base 


RASTER MAP CATEGORY REPORT 


LOCATION: spearfish I Wed Apr 12 09 51 24 1995 


UTM coordinates 

REGION 

north' 2987756.92 

south 2765089 799 

res. 99.98523619 

east' 600898.8994 

west 298182 304 

les. 100 00548246 


REGION CHARACTERISTICS 

REGION 

2987756.92 N 

2765089.80 S 

600898 90 E 

298182.30 W 

ZONE 42 


Longitude Latitude Details 

26.59:58.238147N 

27 00 37.898 189 N 

83.'02.02.996629W 

79.58 58 164221 W 

24.59:23.706679N 

24.59:59 999907N 

82:59:59 412421W 

80:00.00.000003W 


Scale Details 

at northern edge 1 arc second longitude = 99207.680457 meters 
at southern edge 1 arc second longitude = 100911 021934 meters 
at western edge 1 arc second latitude = 110802.149268 meters 
at eastern edge 1 arc second latitude = 110750 608030 meters 
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Table 3.1: Details of Coordinates and Resolutions 
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tion, Reading data from tapes, archiving and retrieving data, documenting procedures 
used to cieate/revise data layer, and maintaining records of data updates 


3.3 Data collection and Map preparation 

3.3.1 Study area 

The maps for Ganga river basin were obtained from the Uttar Pradesh Irrigation 
Depaitment They are the xerox copies of the original maps. These maps are first 
studied caiefully to develop the general understanding of the area. The important 
maps which could be used as a souice for the deiivation of some other related maps 
weie selected for the analysis. 

The area to be studied is the leach of River Ganga between Kanpur and Al- 
lahabad Since it IS desirable that the bouiidaiies of the area are well defined m 
teims of the hydraulic behaviour, it is proposed to have as the upper(Northerii) and 
lower(Southein) boundaiies major canals on either side of River Ganges, viz. ,Purwa 
blanch and Allahabad branch m the noith, and West Allahabad bianch, Fatehpui 
branch and Sasiir Khadeii river in the south; at the west gradient or head contiol 
boundary defined by the hydrographic stations and at the east by the confluence of 
the livers Ganga and Yamuna(Fig 3.1). In order that the boundary conditions are 
well defined in the model, data are required for some area around also The study 
area and the region of data requirements are shown in Fig 3 1 

The area is located between 25 15 and 25 45 N Latitude and 80 and 82 E Longi- 
tude. It is irrigated by the Lower Ganges Canal(LGC) and the Ram Ganga Pariyojana 
in the Ganga Yamuna Doab and fairly extensively in the area to the north of River 
Ganga by Sarda Sahayak Pariyojana. 


3.3.2 Rainfall 


The normal rainfall in the area varies from a maximum ol lOOOrnm near Allahabad 
to around 800mm at Unnao. The coefficient of variation is around SOoccurs in the 
monsoon season from the middle of June to .September with the maximum in .Inly- 
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Raingauge Stations 
Kanpur 
Chakeri 
Unnao 
Haaanganj 
Purwa 
Ghatampur 
Khajuha 
Fatehpur 
Safipur 
Rae-bareli 
Khaga 
Salon 
Sirathu 
Manjanpur 
Kunda 
Bamrauli 
Allahabad 
Phulpur 
Pratapgarh 


Table 3.2' Rain Gauging stations in the Study Region 


August. 

The details of rainfall stations and the monthly rainfall data for the period of 1979 
to 1982 are collected from the U P Irrigation Board and are used for the database 
creation and subsequent analysis. 

The list of rcviiifall stations in tin? study area is pieseiited in 'Bible 3 2 . 
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3.3.3 Ground Water 

Water table levels constitute one of the most important basic data sets for the creation 
and analysis of a database. Adequate ground water level data concerning shallow 
water tabels(piemonsoon and postmonsoon) are collected for different districts in the 
study area. 


3.4 Desk top Scanning 

3.4.1 Introduction 

Over the past 20 years, a great deal of attention has focused on development of tech- 
niques to alleviate the labor intensive process of manual cartographic digitizing as 
tradition ally accomplished using tables and cursors Advances in digital equipment 
foi data entry, in computer processing speeds and in memoiy addressing have in- 
creased the options previously available in digitization. Although manual tables are 
still most widely used, the raster scanner is rapidly becoming a viable alternative A 
scanner takes a very diffeient appioach to capture of digital data from analog origi- 
nals than does the digitizing table. In raster scanning, an oiiginal is captured in its 
entirety in a few minutes or less. A map becomes a massive matrix of reflectivity 
values, each a generalization of a small portion of the original. 

Scanners reduce the amount of time required to produce computer readable file 
and are economically advantageous in systems requiring at least moderate levels ol 
data entry. Throughput time savings from scanning are significantly affected by the 
amount of map preprocessing, and data postprocessing (format conversion [vector- 
ization] and feature labeling) required to create a use able data file. 

fn the recent times, scanning, like many electronic technologies, has declined 
greatly in piire, but laige format pioduction scanning is still (’xpcnsivc and i oin- 
plex enough that manual digitization leitiains the pievaleiit digitizing technique. 
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3.4.2 Raster scanning of cartographic products : 

The capabilities of a scannei are described through specifications of spatial and 
radiometric resolutions. Spatial resolution, refers to size of individual pictuie ele- 
mcnts(pixels) and may be defined by either lows and columns of an output matrix or 
by numbei of ’’dots pei inch"(dpi).either measurement can be converted by the scale 
of the map oiiginal to ground resolution For example if scanning a 1:24,000 scale 
map at 200 dpi the pixel size is 10 feet 

Radiometric resolution describes the range of values the device can discriminate 
for each pixel Radiometric resolutions typically range from two levels of bright* 
ness(bmaiy) to 256 levels of brightness(.255 grey level). 

3.4.3 Desk top Scanning 

’’Desk top publishing’s rapid ascention from an eccentric endeavor to a bonafide soft- 
ware application has earned scanning equipment to the (ore” (O’Alelly, 1988, p 208) 
Today’s fiat bed scanners are alTordable to a wide audience and offer an impiessive 
airay ol hardwaie and software features including: 

• Usei selectable binary and grey /color radiometric resolution(up to 256 levels) 

• User selectable spatial resolution more than 600 dpi. 

• Interactive image sub setting to allow the user to scan a selected portion of the 
oiiginal. 

• Scanning of the original documents up to 8.5 by 14.0 inches and larger. 

3.4.4 About the scanner and output files 

A Hewlett-Packard 9000a ScanJet desk top scanner with associated ’’Scanning Gallery 
Plus 5.0” software was used to scan the photocopies of original maps. 

The software provided with the scanner supports digitization ot documents in a 
two phase procedure, A preview scan of the entire surface is presented in same form 
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on tho computer screen. The user then subsets the image from the screen preview, 
and the final scan recovers data only fiom that poition of the scanning surface, A 
wide variety of options is available for the final scan, and a Tag Image File [TIFF) is 
wiitten to a data storage device. 

Though theie is no single standard format, tag image file contains header followed 
by pixel data. Values in the header reveal number of lovvs and columns in the image, 
and the data aie in different format for binary and for grey scale files As with many 
of the 1 aster formats TIFF files can be extremely large A 256 grey scale radiometric 
file containing 672 rows and 1167 columns requires approximately 784KB of storage 
if data aie packed at 8 bits per pixel If translated into ASCII format { img format) 
the file lequiies approximately 3186KB of meinoiy, more than four times that of a 
TIFF file. However, it is easiest to extiact data from the TIFF file directly 


3.4.5 Spatial and Radiometric resolutions 

Scanning a map at lower resolution (dpi) lesults in significantly shoiter processing 
times as many post processing algorithms progressively thin the lines through peeling 
away edges until lines of a single pixel width result. The opposing view is that a 
certain minimum dpi is requiied to maintain the connectivity of line work and to 
picserve fine detail Trials deteimmed that 150 dpi is the optimal spatial resolution 
to define all the concerned regions on the maps, and that 256 grey scale is much more 
effective than the binary scale for resolving variations among the map colors. Hence 
a spatial resolution of 150 dpi and a radiometric lesolution of 256 grey level has been 
adopted for scanning different maps. 


3.5 Digitization of a Vector map 

3.5.1 Input of spatial data 

Most often vector maps defining the linear features of the region are required for data 
analysis and modeling. The process in which these maps are generated is known as 
digitization. Duiing digitization all featuies are recorded as a series of X and Y- 
coordinates, using a digitizer. In present study the screen itself is a digitizing table 
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and the GRASS bmaiy raster layer seives as the map The GRASS program v. digit, 
a highly interactive digitization and map development program is used for vector 
digitizing, editing, and labeling. The raster map layer is displayed on screen as a back 
drop cell map using the customization mode of the v digit software and using the 
menu based features of the piogiam , details of the linear features of the actual map 
are digitized by careful movement and appropriate clicking of the mouse 

Editing : 

Since digitizing is a tedious opeiation, mistakes will inevitably occui. Editing facilities 
are theiefoie necessaiy, Common eirois occurring during the process of map digitizing 
include’ 

• Failiiie to snap lines together at nodes 

• Ovei- and under shooting of lines 

• Omission of lines oi points. 

• Double recording of lines or points 

• Incoirect feature coding 

• Incoirect location of featmes 

Attribute entry ; 

Feature labels are enteied via keyboard within a spatial menu area reserved, imme- 
diately after a line or an area or a site is digitized. 

3.6 Input of Non-spatial data 

Input of non-apatial data for a map in GRASS is done by editing the coi responding 
vector digital file in dig-att diiectory And after editing this file it is linked with the 
spatial data. This is done by using the utility v, sup port. 
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3.7 Vector to raster conversion 

Since GRASS is a laster based geographical information system the original digitized 
vector data is conveited befoie the analysis is performed The process of vector to 
raster conversion involves the following steps 

The piogram v.to.rast is a conversion utility, which converts vector data to raster 
data 


3.8 Data compression 

Since many of the cells contain the same value as neigh bouring cells, there is consid- 
erable redundancy. Hence significant leductions in size of the data files is achieved 
by using the data compiession option of the program r support 

The list of maps which are digitized and converted into raster format is shown in 
the Table 3 3. 
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S No. 

Map Name 

Title 

1. 

ka.tsp 

Study Region 

2 

ka tsp poly 

Theissen Polygons 

3. 

r sites 

Raingage Stations 

4 

g sites 

Ground water stations 

5. 

pie.gw 

Premonsoon groundwater table levels 

6 

post gw 

Postmonsoon groundwater table levels 

7 

r m 

Monsoon rainfall values 

8. 

i.iim 

Non monsoon rainfall values 

9. 

r t 

Normal rainfall values 

10 

clist 

Districts in the region 

11. 

ganga 

River Ganga 

12 

canals 

DilFeient canals 

13 

cancoin 

Canal command area 

14 

inch a 

Map of India 

15 

up 

Map of Uttar Pradesh 

16. 

blocks 

Irrigation Blocks 

17. 

bl.tia 

Total Irrigated Area( Block wise) 

18 

bl.l 

Net Irrigated Area( block wise) 

19. 

bl.2 

Canal Irrigated Area(blockwise) 

20. 

bl3 

Paddy Irrigated area{blockwise) 

21. 

bl.4 

Sugarcane Irrigated Area(blockwise) 

22. 

bl.5 

Wheat Irrigated Area(blockwise) 


Table 3.3- List of maps pieparecl through v.lo vast 
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Data Analysis and Map Display 


4»1 Introduction 

The most important characteiistic of Geographical Information System is the provi- 
sion of capabilities for spatial analysis functions These functions use the spatial and 
non-spatial attributes in the database to answer questions about the real world 

The database in a G.l S. can be iiserl to simulate certain aspects of reality A 
model may be represented in woids, in mathematical equations or as a set of spatial 
relations displayed as a map 


4.2 Analysis functions 

The objective of data analysis is to extract or query useful information to satisfy the 
requirements or objectives of decision makers at all levels of detail. An important use 
of the analysis is the possibility of predicting what will happen in another location 
or at another point in time This ability provides the possibility to select the best 
possible altci native design. 

The range of analysis procedures have been subdivided into following categories. 
• Retiieval 
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• Reclassification 

• Measurement 

• Overlay 

• Search functions 

• Topographic functions 

• Interpolation functions 

• Other functions 


4.2.1 Retrieval operations 

Retrieval operations on the spatial and non spatial data involve the selective search 
and manipulation, and output of data without the need to modify the geogiaphical 
location of features or to create new special entities. 

Retiieval operations include. 

• Retrieval ol data using geometric coordinates 

• Retrieval of data using symbolic specifications 

• Retrieval of data using conditional and logical statements. 

Theie are numerous retrieval utilities available m GRASS and almost all the util- 
ities are used to in the analysis of maps and generation of reports. Some of these 
utilities are totally interactive(needing a mouse) and some are non-interactive The 
list of these utilities and their brief details are tabulated in Table 4.1. 

For example the utility d profile shows the profile of a particiildr variable along the 
u&ei selected section on the map The utility d lii.-^lo(ji am Aiows the bar/pie (hails ot 
the category values. 

For (Example, monsoon rainfall surface (interpolated from tlic piogram r.siirf idw 
and named r m.surf) and four sections on this map layer are used as Inputs to the 
utility d profile, the results showing the lainfall along the foui sections are shown in 
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program 

Discnption 

d geodesic 

Displays a geodesic line, tracing the shoitest 

distance between two geographic points along a great circle, 

in a longitude/latitude data set 

d histogram 

Displays a histogram in the form of a pie or 
bai chart for a user-specified raster file 

d. legend 

Displays a legend foi a raster map layer in the 
active frame on the graphics monitor 

d ihumblme 

Displays the rhumbline joining two usex- 

specified points, in the active frame on the user’s graphics 

monitor 

d what rast 

Allows the user to interactively query the 

category contents of multiple raster map layers at user- 

specified locations within the current geogiaphic region 

d where 

Identifies the geographic coordinates associated 
with point locations m the active frame on the graphics 

monitor 

d rast.num 

Overlays cell category values on a taster map 
layer displayed to the graphics monitor 

r.cats 

Prints category values and labels associated with 
user-specified raster map layers 

r coin 

Tabulates the mutual occurrence (coincidence) of 
categoiies for two raster map layers. 
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progiam 

Discription 

r.covar 

Outputs a covanance/correlation matrix for user- 
specified raster map layer(s). 

r.deacnbe 

Prints terse list of category values found in a 
raster map layer. 

r info 

Outputs basic information about a user-specified 
raster map layer 

r report 

Reports statistics for raster map layers 

r stats 

Geneiates area statistics tor raster map layers 

r volume 

Calculates the volume of data ’’clumps”, and 
produces a GRASS sitelists file containing the 
calculated centroids of these clumps. 

s.db.iim 

RIM data base management/query interface for 

GRASS sites data. 

r aveiage 

Finds the average of values m a cover map 

within areas assigned the same category value in a user- 

specified base map. 

d.measure 

Measures the lengths and areas of features drawn 
by the user in the active display frame on the graphics 

monitor. 

d, display 

1 Displays profiles of a user-specified raster map layer. 


Table 4 1: Retrieval Opeiations in CiRAS.S4-l 
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program 

discnption 

r buffer 

Creates a raster map layer showing buffer zones 


surrounding cells that contain non-zero category values 

r Cl OSS 

Creates a cross product of the category values 


fiom multiple raster map layers 

r reclass 

Creates a new map layer whose category values 


aie based upon the user’s reclassification of categories in 


an existing raster map layer 

r combine 

Allows category values from several raster map layers to be combined 


Table 4.2 Reclassification Utilities in GRASS4.1 

the Fig 4 10. The block wise irrigated areas are given m map layer hi tia The utility 
d histogiarn is used to present the iriigated areas in each block cis a bai chart (Fig 
4 12) 


4.2,2 Reclassification procedures 

Reclassification procedures involve operations that reassign thematic values to the 
categories of an existing map as a function of the initial value, the position, size or 
shape of the spatial configuration associated with each category, for instance a rainfall 
surface map, leclassified into a map of irrigation blocks involves the process of looking 
at the attribute for a single data layer and assigning an additional attribute, the new 
clciss name. 

Reclassification can also be performed in multiple data layers as part of an overlay 
operation. For example, a sugar cane irrigated area, a paddy irrigated area and a 
wheat irrigated area map can be presented as different data layers. 

The reclassification utilities picsent tor the ruialy.sis in GR/\.,S.S l I and llicMi biud, 
is tabulated in Table 4 2 

r cross and r combine are the reclassification utilities which were lused in the 
present work. The utility r.combine helps us to combine different maplayeis and 
also to overlay them The piogram r.cross generates the cross pioduct of the values 
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Average Monsoon rainfall (as given by program) 

Station 

Rainfall(mm) 

Kanpur 

740.4337156575 

Chakeri 

755 9480371859 

Unnao 

737 1166480962 

Hasanganj 

748 8391180073 

Purwa 

765 0584144645 

Ghatampui 

752 4656925032 

Khajulia 

783.5278414453 

Fatehpur 

787.972397369 

Safipur 

810.6612679858 

Rae-bareli 

822.0424881703 

Khaga 

828.8577208558 

Salon 

845.5810441954 

Sirathu 

846.1989276627 

Manjanpur 

817 8595643596 

Kunda 

853.5689632589 

Bamrauli 

853.2688271784 

Allahabad 

845.6907494936 

Phulpur 

865.8317571776 

Pratapgarh 

868.8675237563 


Table 4.3: Average mon.soon rainfall in difFerent zoiies(Repoit from r.ooss) 
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program 

discription 

d geodesic 

Displays a geodesic line, tracing the shortest 

distance between two geographic points along a great circle, 

m a longitude/latitude data set 

d. histogram 

Displays a histogram in the form of a pie or 
bar chart for a user-specified raster file 

r. report 

Reports statistics for raster map layers 

r. volume 

Calculates the volume of data ’’clumps”, and 
(optionally) produces a GRASS site-lists file containing the 
calculated centroids of these clumps 

r average 

Finds the average of values m a cover map 

within areas assigned the same category value in a user- 

specified base map 

d measure 

Measures the lengths and areas of features drawn 
by the user in the active display frame on the graphics 

monitor 

V area 

Measuies the areas of specified regions on a vector map 


Table 4.4 Mecisurement functions m GRASS4.1 
present in the given map layers. 

4.2.3 Measurement functions 

Measurement of spatial data involves the calculation of distances between points, 
lengths of lines, area and perimeter of polygons, and volumes. 

There are number of utilities which perform the functions of measurement. For 
c'cainple the map layer r.m.suif is used as the upper map to calculate tlie aveiage 
monsoon rainfalls over different districts in the study area (map layer dist) using the 
utility r.average. The results are tabulated in Table 4.10, 

The list of medsuiement functions in GRASS4.i and their brief is given in Table 


4.4. 
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4.2.4 Overlay operations 

Ovei laying of maps result m the creation of a new map where the values assigned 
to every location on that map are computed as a function of independent values 
associated with that location on two or more existing maps Overlaying creates a 
new data set containing new polygons formed from the intersection of boundaries of 
the two or more data layers In addition of creating new polygons based on the overlay 
of the multiple layers, these polygons have multiple attributes, i,e., the attributes which 
were given to each separate overlay before the composition occured. 

Arithmetical and logical overlay are also possible in GRASS such as addition, 
subtraction, division, multiplication of each value m a data layer by the value cor- 
responding to the same location in another data layer Logical overlay involves the 
selection of an area where a set of conditions are satisfied 

The typical GRASS overlay utilities and their brief is listed in Table d 5. The 
utility r. combine is a languauge driven overlaying utility and is used to overlay the 
map layers ganga and cancover showing the River Ganga and the canal command 
area respectively on the study area (map layer ka tsp) The report is shown in Table 
4.6. The net irrigated aieas in each block {map layer bl 1) are overlaid on the raster 
map showing the blocks map layer blocks) using the utility r cross The report is 
shown m Tables 4 7 and 4 8 

4.2.5 Search functions 

Search functions constitute one of the most commonly used neighborhood functions 
They determine the value of each target feature according to some characteristic in 
the neighborhood. 

Typical search functions are the functions which determine the total, average, 
maximum, minimum and statistics siuh a.s mean, standard deviation etc 

The typical GRASS search functions, and their brief are shown in Table 4.9. and 
the resultant reports are tabulated in Tabled. 10, 
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program 

discription 

r combine 

Allows category values from several raster map layers to be 
combined 

r mapcalc 

Raster map layer data calculator, performs arithmetic on 
raster map layers 

r cross 

Creates a new map layer whose category values 

are based upon the user’s reclassification of categoiies in 

an existing raster map layer 

r weight 

A language driven raster map overlay program 

r infer 

Outputs a raster map layer whose category values 
represent the application of user-specified criteria (rules 
statements) to other laster map layers’ categoiy values 


Table 4 5' Oveilaying Functions in GRASS4 1 


MAP 

■ gcover in tsp 

OVERLAY 1 

Non Canal Command Area 

OVERLAY -2— 

Ganga/Yamuna River 

1—2 

River Gang a 

1 :i— 

Canal Command Aica 


Table 4.6* Overlay of command areas and rivers 
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Net irrigated area 

in irrigation blocks(Hectares) 

NAWABGANJ 

6898 

HASANGANJ 

15947 

MIYANGANJ 

8109 

SAFIPUR 

10696 

SIKANDARPUR 

9622 

KALYANPUR 

8433 

BIDHNOO 

12407 

BHITARGAON 

14391 

KANPUR 

11718 

PURKARON 

10236 

BICHCHIA 

17086 

ASOHA 

8925 

PURWA 

11501 

BIGHAPUR 

5090 

SARSAUL 

14090 

HILAULI 

7620 

KHIRON 

9366 

SUMERPUR 

9000 

DEOMAI 

10469 

KHAJUHA 

11707 

SARENI 

9878 

LALGANJ 

7672 

DALMAU 

12381 

SATAON 

10709 

TALYANI 

11667 

HASWA 

11661 

BITHURA 

11882 

RAH I 

1.1993 

JAG ATP UR 

12413 

NATHGARH 

6800 

VIJAIPUR 

6954 

DHATA 

9128 


Tabic 4.7; Report from r. cross 
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SARASWAN 

8156 

UNCHAHAR 

12074 

SALON 

14368 

KALAKANKAR 

9347 

KARA 

5266 

RAMPUR 

12744 

BABAGANJ 

10876 

KUNDA 

9957 

SIRATHU 

7183 

BIHAR 

10945 

ARASWAN 

4400 

MANJANPUR 

4558 

NEWADA 

4906 

KAURIHAR 

5564 

CHAIL 

4539 

KAURIHAR 

5564 

HOLAGARH 

7866 

MAU-AIMA 

9194 


Table 4,8' Net Irrigated Area( in Hectares) 






Data Analysis and Map Display 


46 


program 

discripton 

r average 

Finds the average of values in the cover map within areas as 


signed the same category values in base map 

r watersh.ed4.0 

Watershed basin analysis program 

r basins fill 

Generates a raster map layer showing 


watershed subbzisins 

r cost 

calculates the cost to go to one point to another. 

r drain 

shows the drainage lines 

r thin 

Smoothens the raster map layer 

r grow 

Increases the category zone by one pixel 


Table 4 9: Search Functions in GRASS4.1 


MAP: Reclass of dist.l in tsp (m.ave in tsp) 

Unnao 

759.2609763687 

Kanpur 

752.0110575723 

Rae- Bareli 

818 5627412587 

Fatehpur 

805.2529027288 

Allahabad 

841.9527644379 

Piatapgarh 

8.56 22646091 


Table 4,10: Average monsoon rainfall over districts 
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r surf.contour 

r slope aspect 

r cost 

r.traj 

r. volume 

r.los 

d.rast arrow 

r. watershed 


Table 4 IT Topographic Functions m GRASS4.1 

4.2.6 Topographic functions 

The suiface chaiactenstics such as the slope, relief and form of an area are referred to 
as topography 

In the present study these functions are not used, but the list of these utilities is 
tabulated in Table 4 11. 

4.2.7 Interpolation functions 

Interpolation is the procedure of estimating unknown values at un sampled sites 
using known values of existing obseivations at neighborhood locations. Point based 
interpolation is used to estimate the values at predetermined locations using points 
of known locations and values. The output can be an isohyetal map The utilities 
r surf.idw and s.surf.idw are the programs which geneiate the trend surface from the 
available values. The quality of the interpolation results is a function of precision, 
accuracy, number and distribution of points used in the calculation and the manner 
in which the mathematical function models reality. The choice of appropriate model 
13 , therefore crucial in order to secure reasonable results. The interpolation method 
used in the present study is explained below. 

4.2.8 Inverse squared distance weighting: 

One of the most often used methods of interpolating a value of a variable Z at an 
unvislted point x, i.3 to compute an average value from local neigbboihood, In its 
simplest form, for a regularly spaced data along the transect, the moving aveiage for 
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a point X in tlie centre of symmetric window is computed as 


Z{x) = 1/n ( Z{x,) I (4 1) 

In two dimensions, the same formula would apply, with the sitesr, replaced by 
the coordinate vector X,, 

The size of window has a definite effect in the form of smoothed output Nanow 
windows will emphasize short range variations and broad windows will reduce short- 
range variations m favour of longer range effects. 

As we have alieady noted, observations located close together tend to be more 
alike than observations spaced further apart. It is natural to feel that the contribution 
that a given sample point makes to an average interpolated value at an unvisited site 
should be weighted by a function of the distance between that observation and the 
site So we can compute a weighted moving average 


2(a:) = 5 : A.Z(i,) = l (“IS) 

t = l 

where the weightsA, are given by a:,)). 

A requirement is that <^{d) 00 as d — ^ 00 , which is given by commonly used 

reciprocal or negative exponential functions. Perhaps the most common form of^(d) 
IS the inverse squared weighting 


. _ Er,, z(x,)d:/ 

^1=1 “ij 


(4 3) 


where Xj are the points at which the surface is to be interpolated; usually the 
points lie on a regular giid 


The size of the domain not only affects the aveiage value estimated at a point, but 
also the amount of computer time required for interpolation. Usually the size of the 
domain is set to a minimum and a maximum number of data points in an effoit to 
balance computational efficiency against precision. The number of poinhs used may 
vary between 4 and 12 , but is usually in the range of 6 - 8 , particularly if the original 
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r siiif.idw 

s surf idw 

1 contour 



Table 4.12 Interpolation Functions in GRASS4 1 

data lie on a regular grid, but when the data are irregularly distributed, each inter 
poldtion Will be made using a window of different size, shape and orientation. 

The clumping problem with iriegularly spaced data is solved by giving reduced 
weights to those data points that are scieened fiom the interpolation point by other 
points This is done by dimmishing the weight given to the farther point in proportion 
to the cosine of angle between the data points and the interpolated point 

The interpolated values at the grid points are displayed as greyscale i aster maps 
These inteipolated giid cell values can be used as an overlay in a raster data analysis 

In the present study the monsoon nonmonsoon and annual rainfalls ovei the study 
aiea aie interpolated from the rainfall values at the lain gauging stations using the 
two models s surj idw and r surf.idw The interpolation for the two models differ In 
that 6 surf.idw gave steep profiles not representative of the variation of precipitation, 
hence in this case r surf tdw seems prefeiable The result for total annual precipitation 
IS shown in Fig 4 4 

Pie-monsoon and post-monsoon watei table levels for 131 stations die available. 
They weie used to derive the ground water table surface and contours 


4.3 Map Display 

The maps which were developed using the above mentioned procedures have been 
displayed on a screen using several GRASS display utilities The GRASS display 
tools are tabulated in Table 4 13. The list of the developed Raster, Vectoi and Site 
maps IS shown in Table 4.14. The photographs of some of those maps are shown m 
Fig. 4 1 to Fig 4 12, the details of these figures are: 
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0 The laster map of Study Area ( Gaga River Basin between Kanpur and Alla- 
habad is shown in Fig 4.1 

• The raster map of the six districts in the study region are shown in Fig4 2 

• The 1 aster maps of canal system and the study area are overlaid and is shown 
as a single map layer in Fig 4.3 

• The annual rainfall surface is interpolated and this surface is overlaid on the 
study region and is shown in Fig 4.4 

• The map showing the rainfall Theissen polygons is shown in Fig 4 5 

• The vector map of rainfall contours is overlaid on the study area and the resul- 
tant map layer after appropiiate labelling is shown in Fig 4 6. 

• The vector map of post monsoon ground water table level contours (extracted 
fiom the post monsoon ground watei level interpolated surface using the utility 
r contour) is oveilaid on the map lepresenting the study area. This map is 
shown in Fig 4.7 

• 3D view of topographical surface of the study region is shown in Fig 4 8 

• 3D view of Annual rainfall distribution over the study region is shown in Fig 

4.9. 

• The color photo showing the capability of the utility d ptofile is shown in Fig 

4.10. 

• The rainfall Theissen polygons m the study area with their legends aie shown 
in Fig 4.11 to illustrate the capability of the utility d display. 

• The raster map of the total irrigated areas in each block and the bar chart show- 
ing the values of these areas are shown in Fig 4 12 to illustiate the capabilities 
of the utilities d.lmtoyrnmmxd split viewsli 
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Raster 

d display 

d his 

d rgb 
d. profile 

d 3d 

3d. view sh 

slide.show.sh 

Vector 

d. display 

d vect 

d.vect dig 

Sites 

d display 
d.points 

d sites 

d icons 

Text 

1 

1 

d. display 

d.menu 

d text 

d. label 

d. labels 
d paint labels 
d.legend 

show fonts.sh 

Screen 

Graphics 

d display 
d scale 

d.grid 

d colortable 

d. graph 
d.mapgraph 
d geoflesic 
d.rhumbliiie 

d. histogram 
grass logo,.sh 




O M Hii 

7 




Table 4.13 Map display functions in GR,ASS4.L 




Data Analysts and Map Display 


52 


Map Name 

discription 

ka tsp 

Study Region 

ka tsp poly 

Theissen Polygons 

r sites 

Raingage Stations 

g sites 

Grouad water stations 

pre.gw 

Premonsoon groundwater table levels 

post gw 

Postmonsoon groundwater table levels 

r.m 

Monsoon rainfall values 

r.nm 

Non monsoon rainfall values 

r t 

Normal rainfall values 

clist 

Distiicts in the region 

ganga 

River Ganga 

canals 

Different canals 

cancom 

Canal command area 

India 

Map of India 

up 

Map of Uttar Pradesh 

blocks 

Iriigation Blocks 

bl.tia 

Total Irrigated Area(Block wise) 

bl.l 

Net Irrigated Area(blockwise) 

bl.2 

Canal Irrigated Area( block wise) 

bl,3 

Paddy Irrigated area(blockwise) 

bl,4 

Sugarcane Irrigated Area(bIockwi3e) 

bl,5 

Wheat Irrigated Area(blockwise) 

r. sites. surf 

Rainfall Surface(Total) 

r.nm. surf 

Rainfall Surface(Nonmonsoon) 

r.m. surf 

Rainfall Surface(monsoon) 

pre.gw sutf 

pie nion g w table suiface 

post. gw. surf 

post.mon g.w. table surface 
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Map Name 

Discripton 

r. sites cont 

total rainfall contours 

r.nm.cont 

r {.contours (nonmonsoon) 

r m cont 

r.f contours (monsoon) 

pre.gw cont 

g w t. contours (premon) 

post gw. cont 

g.w.t contours (postmon) 

cross, nm.surf 

cross of nonmonsoon r f and polygons 

cross m.suif 

cioss of monsoon r f and polygons 

cross t.surf 

cross of toal r f and polygons 

pregw.surf 

mtei pointed sueface of premonsoon g w.t 

postgw suit 

interpolated suiface of postmonsoon g w t 

cross tia 

cross of blocks and total irrigated area 

cioss.l 

cross of blocks and net irrigated area 

cross. 2 

cross of blocks and canal irrigated area 

cross 3 

cross of blocks and paddy irrigated area 

Cl 033.4 

cross of blocks and sugarcane irrigated area 

cross. 5 

cross of blocks and wheat irrigated area 

cross 6 

cross of blocks and double iingated area 

cross. 7 

cross of blocks and noncaiial irrigated area 

C10SS.8 

cross of blocks and (paddy-f-wheat) double irrigated area 

cross. 9 

cross of blocks and kharif canal irrigated area 

cross 10 

cross of blocks and rabi canal irrigated area 

cross. 11 

cross of blocks and canal irrigated area(wheat only) 

cross. 12 

cross of blocks and non canal irngsted area(wheat only) 

cancover 

overlay of cancom over gaiiga 

caricovei 1 

overlay of cancom over canals 
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Map Name 

Discription 

gcover 

overlay of cancover over ka 

m ave 

Average rainfall in different districts 

nm.ave 

Average nonmonsoon rainfall in diffeient districts 

post ave 

Aveiage postmonsoon g.w t in districts 

pre ave 

Aveiage premonsoon g.w t in districts 

sr m.ave 

Average monsoon rainfall over the entire region 

sr.iim ave 

Aveiage nonmonsoon lainfall over the entire study legioii 

sr pre ave 

Average premonsoon g w.t. over the study region 

sr post ave 

Average post monsooon g.w.t over the entire study region 

sr t ave 

Average total rainfall over the entire study region 

t.ave 

Reclass of districts 

bl.lO 

Rabi canal irrigated area 

bi a 

Canal irrigated area (wheat only) 

bl.l2 

Non canal irrigated area( wheat only) 

bL6 

Double irrigated area ' 

bl7 

Non canal irrigated area 

bl.8 

paddy+wheat double irrigated area 

bl.9 

Kharif canal irrigated area 


Table 4 14: List of Raster Maps in the Data base 
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Vector files available in 

mapset tsp 

1 

dist 1 

pre. gw cont 

2 

clist.poly 

pie.gw cont. tmp 

3 

gs 

r in.cont 

bl 1 

gw 

r. Sites 

bl 1 tmp 

gw. temp 

i sites cont 

bl.2 

md 

r sites.m 

bl.3 

India 

r sites, nm 

bl.4 

indial 

r. sites nm cont 

bl.5 

ka.grid 

r.sites.nm coutl 

blk 

ka m cont 

raingages 

block. aiea 

ka mapcalc 

isites.cont 

block, tia 

ka nm cont 

test 

blocks 

ka sites cont 

tsp. gw 

canal 

ka.tsp 

tsp pregw 

canal 1 

ka.l&p.poly 

tsp2 gw 

cancom 

ka t&p2 

up 

cover 

kapoly 

V. patch 

cover canal 

poly 1 

veg 

cross 

post.gw 

wells 

(l.cont 

post, gw cont 

wells 1 

dill' cont 

post.gw cont. tmp 

wells3 

clisfc 

pre gw 



Site list files available in mapset tsp 

1 

g.sites 

r sites 

up 

2 

gw 

r.sitesl 

wells 

3 

pregw 

raingages 

wells 1 

c. sites 

pregw surf, land 

tgl 

wells2 


3D viewing parameteis files available in lu.ipset tsp 

ka.pregwsuif.3d 

ka.icont.3d 

ka.rgsiirt.3d 
ka.tsp poly 

kapoly,3d 
kapoly. precont.3d 
kapoly. rcont. 3d 
: kapoly.igsnrf 3d 

pie gw.snif 
piegw..siuf.;}<l 

i. bites surf 


Table 4 15’ List of pieparotl maps 
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Figure 4.1: Raster Map of Study Area 


Figure 4,2: Ra.ster Map of Districts ia the study area 
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Figure 4.3: Raster Map of Different canal systems 



Figure 4.4: Raster Map Interpolated Rainfall Surface 
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Figure 4.5: Raster Map Rainfall Tlieisseii Polygons 



Figure 4.6: Raster Map of Rainfall Contours (monsoon) 
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Figure 4.10: Color Plioto allowing the capability of d.projilc 
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Figure 4.11: Color Photo showing the capability of d.display 


Figure 4.12: Color Photo showing the capability of d.hislogram and split, oicw.sh 



Chapter 5 


Ground Water Balance 


5.1 General 

Tire physical and hydiologic c hatcictenhlics of lire study area discussed in the pre- 
vious cliapLeis piovith' the necessaiy barkgroiiiid itifoiiuation for identification and 
estimation of diirprent c oiiipoireuts of ground water balance The area is bourrded 
by livers as well as canals and tiii'ie ate number of small rivers also They could 
be unpoitant from the point of view of iirteiaction with the ground watei system 
However, the inajoi roiitnbiitioii of ground water lecharge is from rainfall over the 
study area, with its time and space variability The recharge from canal and ground 
water iirigatioii systems as well as field percolation from iiiigated areas also consti- 
tute sources of lec liargc to giomid walei The main soui'ces of extraction from giontid 
water is tliiongh pumping by public and private wells, 

5.2 Ground Water Balance Equation 

Tlie giouiul Wcittn’ balajice oqiiatioii over a apecilied timo peiiod iii geiicial terms (’an 
Ijo stated as * ditinge in storage vohuiie = volume of ^ volmncot tion 

After idcMiti tying tlu^ components ot ie(4)aige and extraction , the above' relationship 
can he vviitteii a,s 

/?)• -I- /2c "h /t*/ 'f* /(/ + = //•’ "f Oq d* k'r T 'h ‘^a "h ^ 
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Gionnd Water Balance 


wheip, 

/?, = Natuial n'diaige Irom precipitation 

= Rediaige clue to seepage of rivers, canals, watei courses, ponds etc 
R; = Recharge fioin irrigation arid other activities 
Ig = Inflow to ground water from other sources 
/, = Influent seepage 

Tp = Withcliawl from ground water storage 

Og = Ground water outflow from basin to other basins 

Et = Evapottanspiiation from gionnd water 

E, = EHluent seepage 

Sg = Change in grouiid water storage 

L = Loss through deep percolation 


5.2.1 Lumped Model : 

Spatially lumped model is consirleied with the entiie area as a single unit and aveiage 
values of parariielois of the area aie taken or assumed for the study. 


5.2.2 Study Period : 

Precipitation is the mam source of supplies. Although its distribution varies greatly 
from period to period, yet the variance from year to year is not so maiked and 
precipitation ovei the area umlei study follows a net annual pattern The peiiod oi 
study loi gioimd water balance is therefore taken as one year. 
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5.3 Recharge Components 

5.3.1 Recharge from Rainfall 

The major souice of techaige of ground water is by piecipitation which infiltrates 
tluough the soil (loiizons The amount of recharge depends on vaiious hydiomateri- 
ological, topographic, humidity, peiineability conditions.The recharge due to rainfall 
IS to be estimated gcneially on the basis of detailed water budget analysis which is 
outside the scope of the study. The following approach based on aieview of available 
inloimation is adopted tor the study The Central Giound Water Board has con- 
ducted a number of studies and has concluded that rainfall lechaige m (J.P can be 
e.stiiuatecl as around 15 percent of the annual lamfall, viz 

Barber and Grtcr have analyzed the data for number of wells in non canal irrigated 
dieas and fitted a lelationship. 

R = 0 27 5; (1 6 Bn- 0.12) 

Where Sy is the specific yield and P,n i!5 the monsoon season rainfall in meters 
The value of Sy ranges fiom 0,12 to 0 15 A value of 0,12 is considered in the present 
analysis. 

For obtaining the value of average monsoon and non monsoon rainfalls over the 
entile region the piograin r average is used. 


5.3.2 Recharge from Canals 

Seepage losses from unlured canals very often constitute about 20 percent ot the 
total diversions. U P. Irrigation Department has earned out the recliaige studies on 
dilFerent canals and these results have been used in the present study to account for 
the canal lecharge component In the water balance equation. 

The details of recharge (district wise) are collected from U.P.Irngation dopaitment 
a,ncl aio implemented in the data base. The average of these Lecharge, is tal,'nlateLl 
from the utility r. average and r.mapcalc. 
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5.3.3 Recharge from field percolation 

The percolation from irrigated fields could not be estimated accurately for different 
ciopa in view of limitation of data and time Generally percolation loss is expected 
to be the highest for paddy and the lowest for wheat and is likely to vary between 20 
to 40%. A value of 20% is adopted in this study. 


5.4 Abstraction Components 

7’he abstiaction from the ground water system mainly consist of the exti action by 
wells Besides this the utilization of ground water directly by free standing trees, and 
consumptive use by vegetation either beneficial or iion-beneficial from shallow water 
table areas is also important For estimation of these components the available data 
for the year 1991 were implemented in the data base and the total recharge has been 
found out using the utility r mapcalc. 


5.4.1 Extraction by Wells 

U P Irrigation Department has provided the data of extractioins from various public 
tube wells, private tube wells, pump sets etc and their estimated draft rates for 
nronisoon and iioiimonsoon seasoirs. The withdraw! values for the year 1991 were 
i near por alec! in this study and are used as an important abstraction 


5.4.2 Evapotranspiration 

The evapotranspiration from the ground water is essentially at potential rate from 
lakes and other water bodies and shallow water table areas near rivers and canals at 
a depth of less than Im below ground level. 

Fiom the water table surface and topography a depth to water table map is 
prepared and the area contributing to potential evapotranspiration is evaluated. The 
total evapotranspiration from shallow water table is determined for the two seasons 
Blaney-Griddh formula is used to calculate the evapotranspiration values for different 
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consumptive use coefficient for different months 



Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 


s 

HU 

0 85 

1.00 

1 15 

1.30 

1.25 

1 10 

0 90 







m 




■■ 






Sugarcane 

0.75 

0 80 

0.85 

0 85 

0.90 

0 95 

1 00 

l.OO 

0 95 

0 90 

0.85 

0 75 


Table 5,1; Consumptive use coefficient k for use in Blaney-Criddle's formula 


crops. 


Eliip 


(5 1 ) 


wheie, u=moiitlily consumptive use 
k— consumptive use factor 

p=:month]y percentage of daytime hours of the year 
t=mean monthly temparature 

Tables. 1. gives the monthly consumptive use crop coefficient k for use in Blaney- 
Criddle formula (after Dastane, 1972) 


5.4.3 Effluent Seepage: 

The ground water is discharged from the basin by seepage to streams. The river could 
be effluent and influent depending on the relative levels of the water table and river 
stages and the intersection of the aquifer by the stream. The water table and liver 
stages fluctuate. The river stage may rise above the watei table during the monsoon 
but would generally be below the water table duiing the uonmonsoon period. As such 
water may be discharged from an aquifer to the stream during nonmonsoon period, 
while for part of period it may seep from the river to aquifer during monsoon period. 
The effluent or influent character of the river therefore varies from season to season 
and from reach to reach. In several rivers the dry weather flow is due to the base flow 
from the giound water and this is also true for rivers Ganga, Yamuna, and Hindou 
etc. In the study area the real situation in a rivei is very complicated. Hence the 
estimation of regenerated base flow except on the basis of actual observations is very 
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clifricult or has to be done with utmost care. In the present study due to the limitation 
of data and time, tlie base flow has been approximated for the ground water balance 
study. 


5.4.4 Soil Moisture storage : 


1 lie watei stored in the unsaturated zone of the soil varies as a result of evaporation 
and tianspiration. This soil moisture conlent can be found out by different methods 
such as giavimetric method , Isometric method etc., But normally the soil moisture 
content is not likely to change form year to year The change in soil moisture content 
is theiefoie neglected in the study 


5.5 Ground Water Balance For Monsoon and Non- 
monsoon Seasons 

A numbei of components of the ground water balance for monsoon and non monsoon 
seasons have been evaluated and are shown in Table 5.2. However influent seepage, 
Ground water inflow or outflow soil moisture storage change were not available for 
the study and so the ground water balance could not be verified. 
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Ciroujjcl Wain components for Monsoon season Heclare-m 


Rcchaigc due to monsoon rainfall 359722 

Recharge due to canals, tubewells, irrigation 390166 


Total Monsoon Recharge 749888 


Monsoon season withdraw) 444785 


Gi ound Water components for Nonmonsoon season Hectare-m 
Available Groundwater Recharge 305013 

Recharge, due to Nonmonsoon rainfall 23251 

Recharge due to canals, tubewells, irrigation 38033 


Total Nonmonsoon Recharge 367187 


Ground Water Diaft 205400 

Evapotranspiration 204648 

Effluent Seepage 80000 


Nonmonsoon withdrawl 490048 


Monsoon seasonal difference 305103* 

Nonmonsoon seasonal difference —123861* 

Ground water storage change 463607* 


Table 5.2; Ground water balance 


’* Diffeieiice between these terms may be due to terms not calculated. 








Chapter 6 


Conclusions 


6.1 Summary 

TIk' Ceographical Inlotmafclon Systems are a poweiful set of tools foi collecting, stoi- 
uig, letileving at will, tianfoim ing and displaying spatial data from the teal world 
Foi a pailicnlai set of pin poses. Additional components of GRASS4.1, a taster based, 
already implemontccl (paitly) GIS were implemented in SPARC LX woikstation at 
Water Resources engineering and management centre, Depaitment of Civil Engineer- 
ing, I I.T., Kanpui. The soitware XGRASS4.1 is tested for its application. Then 
GRASS4 1 and XGRASS4 1 were used foi the hydrologic database creation and anal- 
ysis ol gionnd walei balance components in the ganga river basin between Kanpui 
and Allahalnid 


6.2 Conclusions 

GRASS is a versatile software foi development of spatial database and its analysis. 
In the present .study almost all the utilities of this software have been used for the 
creation, analysis and lepiesentatioii of the database. The created database ran be 
used for further analysis veiy easily.Siiice GRASS has been written in C language with 
facility for shell scripting and implemented in a window environment in a netwoik 
accessible form, it provides an extremely versatile aiul poweiful enviionmeiit foi G [ S. 
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and Hpaiiftl dcTiaion suppoit ^yatems in water resources and Civil Engineering systems, 
jilaiiiiiiig, design, development and management. 

• The Ground water balance study can be used for planning the irrigation devel- 
opmental alternatives 

• The created watei resources data base can be utilized for other hydrologic anal- 
ysis with limited effort, 

• Acemate measuiement of influent and effluent seepage can be done to arrive at 
an exact ground water balance study. 

• Ground watei balance study using a lumped model is a crude method Dis- 
tiibuted model study will give bettei ground water balance in the study area 

6.3 Suggestions for further study 

Piesently papci maps piepared by Irrigation Department aie used as input to G I S 
It’s ju'oduction and collection takes a considerable time These difficulties can be 
avoided by using digital cartographic data or remotely sensed images Moreover 
remotely sensed images tend to decrease the time required for the analysis 

The soil chaiacteristics at different irrigation blocks can be collected to get the 
detailed iiifoiination about the specific yield at those areas and these specific yield 
values can be adopted to get the exact ground water balance for the study region. 

The run-off data may be collected or can be modelled and can be added to the 
present data base to be used for the complete water budget analysis of the study 
legion. 
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